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Summary 

Brassica napus is one of the most important oil crops in the world, and stem rot caused by the 
fungus Sclerotinia sclerotiorum results in major losses in yield and quality. To elucidate resistance 
genes and pathogenesis-related genes, genome-wide association analysis of 347 accessions was 
performed using the Illumina 60K Brassica SNP (single nucleotide polymorphism) array. In 
addition, the detached stem inoculation assay was used to select five highly resistant (R) and 
susceptible (S) B. napus lines, 48 h postinoculation with S. sclerotiorum for transcriptome 
sequencing. We identified 17 significant associations for stem resistance on chromosomes A8 
and C6, five of which were on A8 and 12 on C6. The SNPs identified on A8 were located in a 
409-kb haplotype block, and those on C6 were consistent with previous QTL mapping efforts. 
Transcriptome analysis suggested that S. sclerotiorum infection activates the immune system, 
sulphur metabolism, especially glutathione (GSH) and glucosinolates in both R and S genotypes. 
Genes found to be specific to the R genotype related to the jasmonic acid pathway, lignin 
biosynthesis, defence response, signal transduction and encoding transcription factors. Twenty- 
four genes were identified in both the SNP-trait association and transcriptome sequencing 
analyses, including a tau class glutathione S-transferase (GSTU) gene cluster. This study provides 
useful insight into the molecular mechanisms underlying the plant's response to S. sclerotiorum. 
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Introduction 


Sclerotinia sclerotiorum (Lib.) de Bary is a necrotrophic pathogen 
that has no specific host and infects more than 400 different 
plant species, including sunflower (Hellianthus annuus), chickpea 
(Cicer arietinum), soybean (Glycine max) and rapeseed (Brassica 
napus) (Bolton et al., 2006). Brassica napus is one of the most 
economically important oil crops, and stem rot, caused by 
S. sclerotiorum, results in major losses in yield and quality (Purdy, 
1979; Sharma et al., 2015). 

Ascospores of S. sclerotiorum infect plants through senescing, 
old or injured plant tissue, lower stems and flower petals (Garg 
et al., 2010b). The mycelium then spreads throughout the stem 
and leaves Jamaux et al., 1995) producing enzymes that degrade 
keratin and cell wall components and leading to a build-up in 
oxic substances, such as oxalic acid (Bolton etal., 2006). 
Infection causes pale or dark lesions followed by white patches 
of fungus, which eventually leads to stem rot and the collapse of 
the plant. 

When a plant and pathogen interact, a number of plant genes 
are involved in eliciting a response. Resistance (R) genes detect 
he pathogen attack and activate a signal transduction pathway 
hat in turn triggers the transcription of a large number of 
defence response (DR) genes (Gururani et al., 2012; Singh et al., 
2015). The expression patterns of defence genes in B. napus have 


mainly been studied using B. napus or Arabidopsis microarrays to 
compare resistant versus susceptible B. napus cultivars. Gene 
differential expression levels were found to peak 48 h postinoc- 
ulation, and differentially expressed genes included those 
involved in oxidative burst, molecular transport, signal transduc- 
tion, cell maintenance and development; as well as biosynthesis 
genes for jasmonic acid, ethylene, auxin, glucosinolate and 
phenylpropanoid metabolism (Yang et al., 2007; Zhao et al., 
2007, 2009). In addition, Garg etal. (2013) using protein 
profiling identified a number of B. napus proteins related to 
resistance to S. sclerotiorum, such as pathogenesis-related pro- 
teins and proteins involved in primary metabolism and ethylene 
biosynthesis. However, to date no study has investigated the 
defence mechanism of B. napus using an RNA sequencing (RNA- 
Seq) approach, which is more sensitive and accurate, offers 
higher throughput and is more cost-effective than microarray 
analysis. RNA-Seq has emerged as a powerful tool to examine 
differences in global gene expression. 

Resistance to S. sclerotiorum showed quantitative inheritance 
with additive effects and medium heritability (Fusari et al., 2012). 
Early studies looking for resistance genes in B. napus were 
conducted through QTL mapping, and some pathogen-resistant 
QTLs were identified on chromosomes A1, A2, A3, A5, A6, A8, 
A9, C2, C4, C6, C7, C8 and C9 (Behla et al., 2009; Wei et al., 
2014; Wu et al., 2013; Yin et al., 2010; Zhao and Meng, 2003; 
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Zhao et al., 2006). However, due to the lack of common markers, 
overlapping QTLs were difficult to identify. With the recent 
release of the B. napus genome sequence (Chalhoub, 2014), the 
physical positions occupied by QTLs can be identified by aligning 
the QTL primers with the B. napus genome. Li etal. (2015) 
integrated and compared QTLs for resistance to S. sclerotiorum 
from previous mapping efforts, and conserved QTLs were 
identified on chromosomes AQ (22.5-27.5 Mb) and C6 (29.5— 
36.1 Mb). Nevertheless, QTL mapping based on a biparental 
population is limited due to fewer recombination events. 

With the development of high-throughput sequencing tech- 
nology, genome-wide association study (GWAS) based on linkage 
disequilibrium (LD) has emerged as the main tool for gene 
mapping in plants. GWAS can take advantage of the phenotypic 
variation and historical recombination in natural populations 
without the need to construct a mapping population (Nordborg 
and Weigel, 2008). In recent years, GWAS has been widely used 
to identify genes involved in complex traits in B. napus. Cai et al. 
(2014) identified 43 loci related to yield traits using 451 short 
sequence repeats (SSR) and 740 amplified fragment length 
polymorphism (AFLP) markers in192 B. napus inbred lines using 
the Q + K model. Li et a/. (2014a) found two loci associated with 
erucic acid content, four SNP clusters associated with glucosino- 
late content and two loci associated with seed weight using the 
Illumina 60 K Brassica SNP microarray and 472 B. napus acces- 
sions through GWAS. Lu et al. (2014) identified 26 genes related 
to glucosinolate metabolism using a GWAS method termed 
‘associative transcriptomics’ analysis. In addition, GWAS has been 
used to identify resistance genes in crops infected with stem rust 
(Yu et al., 2014; Zhang et a/., 2014a), tan spot in wheat (Triticum 
aestivum; (Kollers et al., 2014), southern leaf blight (Kump et al., 
2011), head smut in maize (Zea mays; (Wang et al., 2012a), blast 
disease in rice (Oryza sativa; (Wang et al., 2014a) and stem 
canker in B. napus (Fopa Fomeju et al., 2014). However, there 
are no reports of GWAS being used to study resistance to 
S. sclerotiorum in B. napus. 

In this study, we carried out a GWAS of S. sclerotiorum 
resistance in the stems of 347 B. napus accessions genotyped 
with the 60 K Illumina Brassica SNP microarray and identified 
significant associated SNP loci. In addition, we performed 
transcriptome sequencing of resistant (R) and susceptible (S) 
accessions inoculated with S. sclerotiorum. Defence genes, 
expressed in both accessions, and those specific to the R 
accessions only were identified. Among the genes identified in 
both the GWAS and transcriptome analysis, a tau class glu- 
tathione S-transferase (GSTU) gene cluster was included, which 
was verified by quantitative real-time PCR (qRT-PCR). Our results 
reveal candidate genes involved in stem rot resistance and provide 
an insight into the genetic mechanisms of resistance to S. scle- 
rotiorum in B. napus. 


Results 
Phenotypic variation 


We evaluated stem resistance to S. sclerotiorum in 347 B. napus 
accessions during 2013 and 2014, with two replications per- 
formed each year. Extensive phenotypic variations were observed 
in relative susceptibility (Table 1). Relative susceptibility was 
calculated according to the formula: the lesion length on the 
tested accession divided by the control, ‘Zhongshuang 9’. The 
relative susceptibility, stem resistance assessment (RA), ranged 
from 0.58 to 1.36, with an average of 0.918, in 2013 and from 


Table 1 Phenotypic variation in Brassica napus 


Trait Mean SD* Min Max Cvi(%) 
2013RA 0.918 0.144 0.58 1.36 15.73: 
2014RA 0.917 0.144 0.59 1.44 15.71 


*Standard deviation. 
*Coefficient of variation. 


0.59 to 1.44, with an average of 0.917, in 2014. ANOVA of stem 
RA revealed significant differences among the genotypes, but no 
significant differences between the data collected in the two 
years (F = 0.02, P = 0.88) (Table $3). Stem resistance assessment 
in B. napus showed continuous variation and approximated a 
normal distribution (Figure 1), suggesting that resistance was 
controlled by multiple genes. The broad sense heritability of stem 
resistance was 84%, which is consistent with the results of Mei 
et al. (2013) (78%). 


Relative kinship 


SNPs from the 60 K Illumina Brassica SNP array were used as 
queries for BLASTn searches against the B. napus genome; 
10 112 of the 52 157 SNPs were not mapped and were 
excluded. An additional 11 113 SNPs (with call frequencies of 
<90% or minor allele frequencies (MAFs) of <0.05) were 
removed. Therefore, 30 932 polymorphic SNPs were used for 
further analysis (Table S4). 

We evaluated the relative kinships among the 347 accessions 
using SPAGeDi, and found that 80.9% of the kinship coefficients 
between lines were <0.05 and 55% were equal to O (Figure S1). 
Furthermore, 14.6% of the coefficients were between 0.05 and 
0.15. Therefore, the kinship analysis results show that the 
majority of the accessions have a null or weak relationship with 
the other accessions in this population of B. napus. 


Genome-wide association analysis 


Trait-SNP association analysis was performed using the P+K 
model (Price et al., 2006). The QQ plot displayed in Figure 2a 
shows that the model could be used to identify association 
signals. A total of 17 significant associations were detected on A8 
and C6 at P < 3.23 x 107° (P = 1/30932, —logioP = 4.49), with 
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Figure 1 Frequency distribution of stem resistance in Brassica napus. 
2013RA and 2014RA refers to the stem resistance assessment (RA) 
evaluated in 2013 and 2014, respectively. The frequency represents the 
proportion of accessions. 
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Figure 2 Manhatttan and quantile—quantile plots resulting from the 
GWAS results for stem resistance in Brassica napus. (a) Quantile-quantile 
plot for stem resistance; (b) Manhattan plot for stem resistance. The 
dashed horizontal line indicates the Bonferroni-adjusted significance 
threshold (3.23 x 107°). 


Table 2 Markers with significant association for resistance to 
Sclerotinia sclerotiorum 


Chromosome Position SNP Allele MAF P value 

A8 15094063 T/G 0.485 1.32E-05 
A8 15095383 A/G 0.499 2.74E-05 
A8 15096961 C/T 0.477 2.15E-05 
A8 15101546 A/G 0.46 2.34E-05 
A8 15101600 A/G 0.465 2.28E-05 
C3_random 4568807 T/C 0:276 1.83E-05 
C6 30784152 T/G 0.357 2.20E-05 
C6 30791903 A/G 0.388 5.98E-06 
C6 30811671 A/G 0.377 8.93E-06 
C6 30812107 A/C 0.371 9.06E-06 
C6 30812416 T/C 0.38 8.93E-06 
C6 30812429 A/G 0.385 9.06E-06 
C6 30813043 A/G 0.393 9.08E-06 
C6 30813372 A/G 0.356 1.38E-05 
C6 30813377 T/C 0.354 2.30E-05 
C6 30813548 T/C 0.364 1.66E-05 
C6 30815877 A/G 0.38 8.93E-06 
C6 31344734 T/C 0.392 3.11E-05 


a false discovery rate (FDR) of 1.2% and 0.8% (Figure 2b). The 
location and MAF of these 17 loci are listed in Table 2. In 
addition, one SNP was found on the C3-random chromosome, 
but this was not analysed further. 

Haplotype block structures were investigated for chromosomes 
A8 and C6. A total of 31 and 13 haplotype blocks were found on 
A8 and C6, respectively. Significant SNPs on A8 were located in a 
haplotype block of 409 kb (from 15.01 Mb to 15.42 Mb; 
Figure 3). For C6, the region of significant association ranged 
from 30.78 Mb to 31.34 Mb, which was consistent with previous 
QTL mapping results (Figure 4) (Wu et al., 2013; Zhao et al., 
2006). 
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Figure 3 Genomic location of SNP loci associated with stem resistance 
(15.09-15.10 Mb) and pairwise LD estimates in A8 haplotype block 
(15.01-15.42 Mb). The red arrows denote the significantly associated 
SNPs located in the haplotype block. 


Prediction of stem resistance 


The ability to predict stem resistance using genome-wide markers 
was assessed. The predictive ability represented the correlation 
between the predicted and measured field data. We used 20- 
90% of the training population and all of the 30 932 SNP 
markers, and found that the predictive ability and standard 
deviation increased with the size of the training population. 
Higher predictive ability and lower standard deviation were 
optimum; therefore, 60% of the training population was selected 
based on both predictive abilities and standard deviation 
(Table S5). 

When using 60% of the training population (145 accessions), 
and all SNP markers for the prediction, the predictive ability was 
only 0.27, indicating that these makers were not suitable for stem 
resistance prediction. In addition, 300 and 1200 randomly 
selected markers were also used and the predictive ability of 
both sets of markers was low (0.17) (Table S5, Figure 5). Markers 
significantly associated with the trait could be detected in GWAS 
by setting a certain P value, such as P < 0.001, P < 0.01 and 
P< 0.05 (Zhang etal., 2014b). We then selected significant 
association markers identified in GWAS with a certain threshold 
for the prediction, including markers with P < 0.001 (69 SNPs), 
P < 0.01 (288), and P < 0.05 (1203), and the predictive abilities 
were 0.58, 0.77 and 0.80, respectively. Therefore, loci detected 
using GWAS improved the predictive ability of stem resistance. 


Transcriptome sequencing 


RNA from five R and five S genotypes was pooled 48 h 
postinoculation to generate R4s and Sug, while mock-inoculated 
B. napus was dubbed Ro and So. A total of 22.0, 20.2, 22.1 and 
19.2 million raw sequence reads were generated from four 
libraries Of Ro, Rag, So and Sag (Table S6). After the removal of 
low-quality reads and adaptor sequences, 18.4, 16.9, 18.6 and 
16.0 million clean reads were obtained, respectively. Eighty per 
cent of the reads were successfully mapped to the reference 
genome, B. napus, in which unique reads that only aligned with 
one loci of the genome accounted for about 89%. 

A total of 6,821 significant differentially expressed genes 
(DEGs) were identified in the two genotypes based on the 
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Figure 4 Comparison of the GWAS results with those of previous QTL 
mapping reported in Zhao et al. (2006) and Wu et al. (2013). The markers 
that lie within the confidence interval of the Sclerotinia stem resistance 
QTL are in blue for Zhao et al. (2006) and red for Wu et al. (2013). 
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Figure 5 Predictive abilities with different training populations and SNP 
markers. ‘All’, means all of the 30932 SNP markers; ‘Random 300’, 
randomly selected 300 markers; ‘Random 1200’, randomly selected 1200 
markers; ‘P < 0.001, P < 0.01, P < 0.05’, significant association markers 
identified in GWAS with the certain P value, P < 0.001, P< 0.01 and 

P < 0.05, respectively. 


criterion, |log2(FPKMas/FPKMo)|=2 and FDR < 0.01. Of the 5384 
DEGs in the R genotype, 2356 (43.8%) were up-regulated and 
3028 (56.2%) were down-regulated. A total of 5386 differen- 
tially expressed genes were identified in the S genotype, 41.4% 
(2229) of which were up-regulated and 58.6% (3157) of which 
were down-regulated (Figure S2a). In addition, for the up- 
regulated genes, 1784 were common to both genotypes (R and 
S), 572 genes were specific to the R genotype only and 445 were 
specific to the S genotype. For the down-regulated genes, 2,166 
were commonly down-regulated (both R and S), 862 were 
specific to R and 991 specific to the S genotype (Figure S2b). 
Finally, hierarchical cluster analysis was performed based on the 
expression levels of the DEGs (Figure S2c). The R and S genotypes 
had similar expression patterns, but the main differences in gene 
expression were quantitative. 


Transcription factors (TFs) play an important role in plant 
resistance response (Seo and Choi, 2015). In this study, 495 
(7.3%) genes encoding TFs were identified as being differentially 
expressed in the two genotypes, of which 84 were differentially 
expressed only in the R genotype, 133 only in the S genotype and 
278 in both genotypes. These genes were divided into 36 
families, including ERF (64), MYB (63), NAC (61), bHLH (55), 
WRKY (48) and C2H2 (37) (Tables 3 and S7). Nearly all ERF and 
WRKY genes were up-regulated, while 78% of the bHLH genes 
were down-regulated, and MYB, NAC and C>H2 genes were up- 
and down-regulated to similar extents. As numerous ERF and 
WRKY genes were found to be mostly up-regulated, our results 
indicate that these two families play a role in the response of 
B. napus to S. sclerotiorum. 

Mitogen-activated protein kinase (MAPK) is a type of serine/ 
threonine protein kinase, and the MAPK cascade signalling 
pathway is important for plant defence response (Zhang and 
Klessig, 2001). In our study, genes encoding for the common MAPK 
cascade signalling pathway MAPKKK19-MKK9-MPKS5 (Sun et al., 
2014) were up-regulated in both genotypes, suggesting that this 
pathway is associated with response to S. sclerotiorum in B. napus. 


Functional classification of common and specifically 
expressed genes 


To monitor the gene expression pattern after infection, Gene 
Ontology (GO) enrichment analysis of commonly expressed genes 
in both genotypes (R and S) was conducted (Figure S3). The 
analysis showed that genes involved in metabolism, cellular 
processes and defence response were up- or down-regulated 
after infection with S. sclerotiorum in B. napus. Most genes 
located in organelles, mainly the chloroplast, were down- 
regulated in both genotypes, suggesting that photosynthetic 
damage was occurring during infection. In addition, KEGG 
pathway analysis of genes that were differentially expressed in 
both genotypes (common genes) was performed (Table S8). The 
KEGG analysis showed that after B. napus was infected with 
S. sclerotiorum, photosynthesis, glyoxalic acid and carbon meta- 
bolism were suppressed, while the biosynthesis of secondary 
metabolites such as sulphur metabolism (especially GSH and 
glucosinolates) and 2-oxocarboxylic acid metabolism was acti- 
vated. 

To decipher the resistance mechanism of B. napus, we 
functionally classified and determined the metabolic pathways 
of 1434 genes differentially expressed only in the R genotype 
(Figure S2b) using MapMan. Genes involved in the jasmonic acid 
(JA) pathway, lignin biosynthesis, signal transduction, defence 
response and encoding transcription factors were identified and 
divided into groups. Group one contains specific genes related to 
the JA pathway and lignin biosynthesis (Table S9). 

The second group contains genes involved in signal transduc- 
tion including MAPK cascade signalling pathway, calcium sig- 
nalling and G protein-coupled receptors. The 12 specific genes 
involved in calcium signalling and 5 genes encoding for G protein- 
coupled receptors were all induced. Genes involved in the 
signalling module MEKK1-BnaMKK4-BnaMPK3-BnaWRKY26 
(Liang et al., 2013) were all up-regulated in the R genotype. 
Another gene family, the MLO genes, MLO7 and MLO71, was 
down-regulated. Previous research shows that loss-of-function of 
MLO results in broad-spectrum resistance to powdery mildew in 
Hordeum vulgare (barley) (Piffanelli et a/., 2004). 

The third group includes a total of 15 genes involved in defence 
response, such as chitinase genes (CHI and CHIB), BGLU16, seven 
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Table 3 Differentially expressed genes belonging to transcription Table 3 Continued 
factor families 
TF family Total number Change Common R specific S specific 
TF family Total number Change Common R specific S specific ——————— SSS a 
SBP 5 Up 0 o 0 
AP2 6 Up 2 0 0 Down 2 Fi 1 
Down 0 3 1 SRS 2 Up 0 0 0 
ARF 6 Up (0) 0 0 Down (0) 1 1 
Down 4 0 2 TALE 10 Up 0 (0) 0 
ARR-B 5 Up (0 (0) 0 Down 6 3 1 
Down 1 3 1 TCP 1 Up 0 0 1 
B3 3 Up (0) 0 0 Down (0) (0) 0 
Down o 1 2 Trihelix 7 Up 2 1 3 
BES1 1 Up 1 (0) 0 Down 1 (0) 0 
Down 0 0 o WOX 4 Up 2 0 0 
bHLH 55 Up 4 4 4 Down 2 (0) 0 
Down 20 1 12 WRKY 48 Up 31 3 2 
bZIP 13 Up 1 1 Down 4 1 7 
Down 6 0 4 ZF-HD 2 Up 0 0 0 
C2H2 37 Up 15 1 3 Down 1 1 0 
Down 6 6 6 
cu 7 -P p Q H disease-resistant protein genes and five receptor protein genes 
—— 2 ? (RLP6, 32, 41, 54 and 56). 
CO-like 6 =P 0 9 & Group four contains genes encoding transcription factors. A 
Down á 1 3 total of 82 TFs were found in the R genotype, of which 29 were 
GRP 1 =P 7 9 up-regulated in response to inoculation, including seven ERFs, five 
me ‘ n MYBs (two MYB108s, two MYB37s and MYB84), five NACs, four 
DBB ; = : : : bHLHs and three WRKYs (two WRKY61s and WRKY8) (Table $9). 
own 
Dof 10 Up 3 0 0 Combined analysis of GWAS and transcriptome 
Down 3 1 3 sequencing 
ERF 64 P 37 7 I The candidate gene regions on chromosomes A8 and C6 as 
Down 1 3 determined by GWAS analysis were 15.01-15.42 Mb and 30.78- 
pelike j -P 3 0 31.34 Mb, respectively. When we combined these results with 
Down d 1 A those of transcriptome sequencing analysis, 21 genes that were 
GATA 7 “P % E 2 differentially expressed in the R or S genotype were found in this 
Bowe 2 7 interval. Eight of these genes were on A8 and 13 were on C6, 
ad 1 an ° o including four genes with unknown function (Figure 6, Table 4). 
Down: 1 o o In addition, a GSTU gene cluster, BnaC06g31020D (W22), 
GRAS 8 uP 9 a i BnaC06g31030D (W23) and BnaC06g31040D (W24), that 
Down 2 1 4 occurred 330 kb downstream of the peak SNP on C6 was up- 
HE 18 T : a o regulated about 64-fold. Both BnaC06g31020D and 
Down 10 3 3 BnaC06g31030D showed homology with GSTU12 from Ara- 
HSF 8 vp > 0 0 bidopsis thaliana, and BnaC06g31040D was homologous with 
Down 1 0 Z AtGSTU11 (Table 4, Figure 6). These genes have been shown to 
ED 3 -P 1 : e play an important role in the defence response of plants (Jiang 
Down l a ! and Deyholos, 2009; Liao et al., 2014). 
MIKC 6 Up 1 o 2 
Down 1 0 2 Verification of transcriptome sequencing 
me a =P 19 3 E To validate our transcriptome sequencing results, qRT-PCR was 
ici 4 8 used to determine differential gene expression for the genes that 
MYB_related 10 =P 0 i 7 were commonly identified. The real-time PCR results confirm and 
Down 4 i 3 verify the expression patterns obtained with transcriptome 
NAC si na 23 2 t sequencing (Figure 7). 
Down 13 6 3 
MEND f Up i o p Discussion 
Down 1 0 0 
Nin-like 1 Up 0 0 0 Culture of Sclerotinia sclerotiorum and the screening 
Down 0 o 1 method for stem resistance 
RAY 1 -P 0 9 1 The isolate of S. sclerotiorum used in this study (kindly provided 
Down 0 0 0 


by Dr. Mei) was collected from the field, incubated and 
grown in the laboratory. Different B. napus accessions vary in 
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Figure 6 The differentially expressed genes in 
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their pathogenicity and resistance to different isolates of 
S. sclerotiorum (Garg et al., 2010a). In this study, the resistant 
accessions could be correctly and effectively identified using the 
same criteria as would be expected when there is a single isolate 
infection of S. sclerotiorum. If mixed populations of S. sclero- 
tiorum were used, the efficiency of screening for resistance 
would have been reduced and GWAS mapping could not have 
been performed. The QTL identification for resistance to 
S. sclerotiorum of B. napus also used the isolate of S. sclerotio- 
rum (Mei et al., 2013; Wu et al., 2013). However, there are still 
biological limitations when using a single isolate, in that it 
cannot represent the diversity that would generally be observed 
in the field. For example, a number of S. sclerotiorum isolates 
can be found in a field. Sexton et al. (2006) identified 55 unique 
S. sclerotiorum haplotypes with genotypic diversity ranging from 
0.40 to 0.71 in a rapeseed field. Furthermore, the aggressive- 
ness of the different isolates was variable due to physiological 
specialization (Taylor et al., 2015). Due to this variability, it 
would be difficult to identify all of the QTLs in B. napus that are 
associated with resistance to all of the S. sclerotiorum isolates. 
Wei etal. (2014) showed that QTLs identified using the 
detached stem inoculation assay overlapped with those found 
in a diseased field; however, some QTLs were found in the field 
but not in the laboratory. 

In this study, the detached stem assay was used to assess the 
resistance of B. napus to S. sclerotiorum. Previous methods for 
screening stem resistance to S. sclerotiorum have included petal 
inoculation, toothpick inoculation and plug inoculation in the 
field (Yin et al., 2010; Zhao and Meng, 2003). However, the 


position of differentially expressed genes within 
IIl the interval. Bottom, the expression levels in the 
resistant (R) and susceptible (S) genotype. Red 
denotes up-regulated genes; green represents 
down-regulated genes; and gray represents the 
missing value. 


W22 W23 W24 


efficiency of these methods greatly depends on environmental 
conditions, such as rain and/or temperature. Screening could be 
better controlled in the glasshouse, but this involves using 
expensive infrastructure to maintain temperature and humidity. 
In addition, field and glasshouse experiments for screening are 
fairly time-consuming. The detached stem assay is a flexible, fast 
and reliable identification method that can easily be performed in 
a controlled uniform environment (Mei et al., 2012). Very little 
space was needed, and a large-scale estimation, about 150 
stems, only needed a 2-m? space (Mei et al., 2012). Furthermore, 
there was no correlation between lesion length and the diameter 
of stems (Mei et al., 2012), which meant that stem resistance 
could be accurately evaluated regardless of different sizes. 

The PDA media that S. sclerotiorum is cultured on affects its 
virulence. Guo and Stotz (2007)showed that using PDA relative to 
minimal medium, disease symptoms occurred 1 day earlier in 
leaves of Arabidopsis. This is advantageous when using the 
detached stem inoculation technique, as it allows for quick and 
accurate evaluation of stem resistance with minimized water loss 
from the stem. In addition, using isolates with increased virulence, 
it may be easier to identify the most significant QTLs. As such, we 
identified two QTLs on A8 and C6 in this study. Future genetic 
research should take into consideration the virulence of the 
pathogen used for screening. Minimal media (relative to PDA) 
would be better used for inoculation trials to identify how an 
isolate, with different virulence, performs in plant-pathogen 
interactions. Furthermore, additional QTLs can also be identified. 

In this study, we wounded the plant stem to allow S. sc/ero- 
tiorum infection. The pathogen does not need to penetrate the 


© 2015 Society for Experimental Biology, Association of Applied Biologists and John Wiley & Sons Ltd, Plant Biotechnology Journal, 1-13 


Resistance to Sclerotinia stem rot in B. napus 


Table 4 Common candidate genes identified by GWAS and transcriptome sequencing 


Arabidopsis R logs S log2 
Code Gene homologue Description (fold change) (fold change) 
W1 BnaA08g19770D AT1G26560 Beta glucosidase 40 (BGLU40) —4.86** 
W2 BnaA08g20100D AT1G27110 Tetratricopeptide repeat-like superfamily protein 8.3274% 
w3 BnaA08g19850D AT1G35350 EXS family protein 5.04% ** 4.58** 
W4 BnaA08g19990D AT1G26920 Unknown protein —5.04** —3.88* 
W5 BnaA08g20190D AT1G27210 ARM repeat superfamily protein —3.95** =5.67*"* 
W6 BnaA08g20360D AT1G23560 Domain of unknown function (DUF220) 3.7'1* 4.19** 
W7 BnaA08g19560D AT1G26100 Cytochrome b561 —4.19* 
ws BnaA08g19580D AT1G26120 soprenyl cysteine methylesterase-like 1 (ICME-LIKE1) —5.18* 
w9 BnaC06g30020D AT1G68430 Unknown protein 3.67* 
W10 BnaC06g30560D AT1G69390 Homologue of bacterial MinE 1 (MINE1) —3.88* 
W11 BnaC06g29990D AT1G69160 Unknown protein —3.62* —3.94** 
W12 BnaC06g30030D AT1G68440 Unknown protein 3.58* 3.05* 
W13 BnaC06g30160D AT1G68560 Alpha-xylosidase 1 (XYL1) —4.37** —4.91*** 
W14 BnaC06g30180D AT1G68570 Major facilitator superfamily protein (MFS) —4*** —4.5** 
W15 BnaC06g30250D AT4G39670 Glycolipid transfer proteinfamily protein (GLTP) 6.48*** 5.237** 
W16 BnaC06g30370D AT1G68740 PHO1 —4.04** —4.53** 
W17 BnaC06g30470D AT1G26450 Carbohydrate-binding X8 domain superfamily protein —5.84*** —7.09*** 
W18 BnaC06g30580D AT1G69420 DHHC-type zinc finger family protein —5.39*** —3.13* 
W19 BnaC06g30610D AT1G33590 Leucine-rich repeat (LRR) family protein 6.78*** 701 
W20 BnaC06g30040D AT1G68470 Exostosin family protein —3.62* 
W21 BnaC06g30280D AT1G68590 Ribosomal protein PSRP-3/Ycf65 —4.08** 
W22 BnaC06g31020D AT1G69920 Glutathione S-transferase TAU 12 (GSTU12) 5.96*** 6.05*** 
W23 BnaC06g31030D AT1G69920 Glutathione S-transferase TAU 12 (GSTU12) 8.31*** 8.57*** 
W24 BnaC06g31040D AT1G69930 Glutathione S-transferase TAU 11 (GSTU11) 8.96*** 7.85%*** 


REP 2 0.0015 **P <0.01; *P < 0.05. 


Log, (fold change) 


Figure 7 qRT-PCR validation of the expression 
patterns of 24 genes identified both in GWAS and 
transcriptome sequencing. 


natural structural defence barriers during infection. This provides 
favourable conditions for pathogen attack, but may circumvent 
any mechanical defence mechanisms of the plant or any defence 
mechanisms specific to the plant epidermal tissue during the 
initial infection. As a consequence, some defence-associated QTLs 
might not be detected. Many similar QTLs were detected in 
wounding and nonwounding tests examining S. sclerotiorum 
resistance in common bean, but one specific QTL probably related 
to initial defence was identified in the nonwounding assay (Soule 
et al., 2011). Therefore, future research could also include 
nonwounding assays to identify QTLs associated with the initial 
defence mechanism. 


The identification of candidate genes 


S. sclerotiorum is one of the most damaging agricultural dis- 
eases in China and can lead to considerable yield loss in B. napus. 
Previous studies into the resistance of B. napus to S. sclerotiorum 


wi w2 w3 w4 WS we w7 W8 W9 W10 W11 W12 W13 W14 W15 W16 W17 W18 W19 W20 W21 W22 W23 W24 


were mainly conducted using QTL mapping, but in this study, we 
used GWAS implemented with GenABEL software for this 
purpose. We identified 17 SNPs that had significant associations 
with S. sclerotiorum response on chromosomes A8 and C6. 
Identification of SNPs on C6 (30.7-31.3 Mb) is consistent with 
previous QTL mapping results (Li et al., 2015; Wu et al., 2013; 
Zhao et al., 2006). Zhao et al. (2006) evaluated S. sclerotiorum 
resistance in the petioles of B. napus DH lines, and a QTL was 
repeatedly detected in the region of 23.1-35.4 Mb on C6. Wu 
et al. (2013) detected a QTL from 30.3 Mb to 36 Mb using 
B. napus stems inoculated with mycelial agar plugs. Li et al. 
(2015) also identified conserved stem resistance QTLs on C6 from 
29.6 to 36.1 Mb through integration analysis of resistance QTLs. 

A comparison of our GWAS and transcriptome sequencing 
results revealed 21 common genes and a GSTU gene cluster, 
which were verified by qRT-PCR. Eight common candidate genes 
found on chromosome A8 were analysed: BnaA08g19770D, 
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BnaA08g20100D, BnaA08g19850D, BnaA08g19990D, 
BnaA08g20190D, BnaA08g20360D, BnaA08g19560D and 
BnaA08g19580D. On chromosome C6, a total of 13 candidate 
genes were found (Figure 6, Table 4). 

Two of these 21 genes encoded for glucosidases, 
BnaA08g19770D (W1) which is homologous to At? g26560 
(BGLU40) and BnaC06g30160D (W13) which is homologous to 
XYL1. Zhao et al. (2007) found that BGLU40 was expressed 48 h 
after inoculation with S. sclerotiorum in susceptible B. napus. 
Glucosidases catalyse the cleavage of the glucoside plant defence 
molecules, such as cyanogenic glucosides and glucosinolates, 
during the interaction between plant and pathogen (Morant 
et al., 2008). Gucosidases have also been suggested to play a role 
in activating phytohormones (Lee et a/., 2006). Increased expres- 
sion of BGLU7 has been demonstrated to be associated with 
resistance to the necrotrophic pathogen, Cochliobolus heterostro- 
phus in maize (Simmons et al., 2001). The resistance has been 
shown to not depend on the antimicrobial compound or 
hydroxamic acids, but operates via other mechanisms (Simmons 
et al., 2001). In our study, BGLU40 expression decreased 29.04- 
fold in B. napus which may lead to increased susceptibility in 
plants infected by S. sclerotiorum. 

XYL1 (BnaC06g30160D), which encodes for «-xylosidase, 
has been shown to catalyse the hydrolysis reaction of xylogu- 
can ologosaccharides, a member of glycoside hydrolase (GH) 
family 31 (Henrissat and Bairoch, 1993). Miedes et al. (2014) 
showed that an increase in xyloglucan content increased 
resistance to pathogens in Aradidopsis. In Arabidopsis, the 
xly1-2 mutant shows enhanced resistance to the necrotrophic 
fungus P. cucumerina due to the structural changes in 
xyloglucan (Delgado-Cerezo etal., 2012). However, in our 
study, the expression level of XYL7 in the S genotype 
decreased (30.06-fold) more than that of the R genotype 
(20.68-fold). Therefore, the role of glucosidases in the defence 
response of B. napus to S. sclerotiorum may need to be 
investigated further. 

The common gene BnaC06930470D (W17) found on C6 
encodes for the pathogenesis-related (PR-2 family) gene, B-1,3 
glucanase. PR proteins are strongly induced by pathogens or 
related factors (Stintzi et a/., 1993). B-1,3-glucanase inhibits the 
growth of fungi by hydrolysing the fungal cell wall, which 
contains the primary polysaccharide, B-1,3 glucan. In addition, 
B-1,3-glucanase triggers the release of oligosaccharide elicitors, 
which activate the production of antifungal compounds such as 
phytoalexin and other PR proteins (Keen and Yoshikawa, 1983; 
Klarzynski et al., 2000). B-1,3-glucanase has also been shown to 
be responsible for the decreased infection of S. sclerotiorum in 
B. napus (Fernando et al., 2007). When the B-1,3-glucanase 
gene (bgn13.1) from the fungus Trichoderma virens-10 was 
transformed and overexpressed in B. napus, the transgenic lines 
inhibited hyphal growth of S. sclerotiorum (Kheiri et al., 2014). In 
our study, the expression of BnaC06g30470D in S genotype 
decreased (136.24-fold) more than that in R (55.71-fold), which 
might increase susceptibility to B. napus. 

The gene (BnaC06g30610D, W19) encodes for a leucine-rich 
repeat (LRR) family protein which was induced in both genotypes 
in our study, and this finding was verified by qRT-PCR. The LRRs 
are important conserved domains in R genes and have been 
shown to be involved in the interaction between plants and 
pathogens (Kohler et al., 2008; Porter et al., 2009). 

The transcription factor BnaC069g30580D (W18) that encodes 
for the DHHC-type zinc finger protein (ZFP) was significantly 


down-regulated in the R genotype (41.93-fold, P < 0.01). TFs 
play an important role in plant resistance response (Seo et al., 
2015). ZFPs can be classified into several different types such as 
C2H2, C2C2, CCCH, GATA and DHHC, some of which have been 
found to be involved in abiotic and biotic stress (Cui et al., 2002; 
Jan et al., 2013). There are relatively few DHHC types of ZFPs, 
and very little is known about their function. It has been shown 
that DHHC-type ZFPs were repressed by Magnaporthe oryzae in 
rice (Li et al., 2014b). However, the mechanism for this is still 
unknown. In addition, many of the other TFs identified in this 
study belonged to ERF and WRKY families and most were up- 
regulated upon inoculation with S. sclerotiorum (Tables 3 and 
S7). The ERF and WRKY transcription factor are known to be 
involved in biotic and abiotic stress response (Chen et al., 2013a, 
b; Son et al., 2012; Wang et al., 2014b; Yang et al., 2009; Zhao 
et al., 2007). In Arabidopsis, ERF5 negatively regulates chitin- 
induced innate immunity response (Son et al., 2012). 

Resistance genes have been shown to generally cluster in the 
genome (Graham et al., 2002). In this study, a GSTU gene cluster 
was identified 330 kb away from the peak SNP on chromosome 
C6; BnaC069g31020D and BnaC06931030D are homologues of 
GSTU12 (AT1G69920), and BnaC06g31040D is a homologue of 
GSTU11 (AT1G69930). All of these genes were up-regulated in 
both genotypes upon S. sclerotiorum infection. GSTU plays an 
important role in the response to biotic and abiotic stress 
(including S. sclerotiorum) as cellular protectants, reducing free 
radicals and reactive species (Jiang and Deyholos, 2009; Liao 
et al., 2014). In previous studies, following infection with 
S. sclerotiorum in B. napus, the homologues of AT/G69920 
(GSTU12) and AT1G69930 (GSTU7 1) were also up-regulated and 
GST proteins were expressed in resistant B. napus (Garg et al., 
2013; Zhao et al., 2007, 2009). The GSTU gene cluster identified 
in this study confirms their importance in response to S. sclero- 
tiorum by possibly acting as an antioxidant defence mechanism. 

Another gene, BnaC./GMT5.a, an indole glucosinolate methyl- 
transferase, located on C6 at 35.5 Mb, about 4 Mb away from 
our peak SNP, has previously been identified and cloned as a 
candidate gene for resistance to Sclerotinia (Wu et al., 2013). RT- 
PCR analysis showed that no expression was observed in the 
susceptible B. napus lines, and increased expression was observed 
in the resistant (R) lines 24 to 96 h postinoculation (Wu et al., 
2013). However, in our study, the expression of BnaC.IGMT5.a 
(BnaC06g37610D) was also increased 48 h postinoculation in 
both R and S genotypes. 

It has also been shown that the presence of susceptibility 
factors (genes) in susceptible genotypes could increase suscepti- 
bility to necrotrophic fungal pathogens in wheat (Liu et al., 
2006). Therefore, it could be argued that resistance may be 
governed by the absence of a susceptibility gene in the resistant 
genotypes as opposed to the presence of a resistance gene. 
Specific susceptibility genes in the S genotype were detected in 
this study; however, the presence of these genes needs to be 
validated in future studies. To follow up on the genes identified in 
this study, we aim to clone the above candidate genes in order to 
elucidate the differences observed between the R and S geno- 
types in order to confirm their roles in the response to Sclerotinia. 


Transcriptome sequencing 


Zhao et al. (2009) found that the greatest number of differentially 
expressed genes between the resistant B. napus variety Zhongyou 
821 and susceptible cultivar Westar occurred 48 h after inoculation 
with S. sclerotiorum. Therefore, we conducted our transcriptome 
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sequencing analysis on B. napus stems 48 h after inoculation. 
Genes involved in photosynthesis, glyoxalic acid and carbon 
metabolism pathways were mostly down-regulated upon inocula- 
tion. Soybean and tomato had a similar response to Pseudomonas 
syringae; genes involved in photosynthesis were down-regulated 
after inoculation (Berger et al., 2004; Zou et al., 2005). 

Conversely, we found that the biosynthesis of secondary 
metabolites, sulphur metabolism (especially of GSH and glu- 
cosinolates) and 2-oxocarboxylic acid metabolism were 
increased, suggesting that the plant immune system was 
activated and plays a key role in the defence response to 
S. sclerotiorum. A total of 81 up-regulated genes were involved 
in sulphur metabolism, and sulphur is an essential macronutri- 
ent for plant growth and development. Recent research into 
sulphur-induced resistance (SIR) or sulphur-enhanced defence 
(SED) (Rausch and Wachter, 2005) shows that GSH and 
glucosinolates induce SIR/SED. GSH is a major antioxidant 
molecule in plants and is reported to be involved in the 
interaction between plants and microbes, including necro- 
trophic pathogens (Ball et a/l., 2004; Frendo et al., 2013; Parisy 
et al., 2007). GSH has been shown to increase after inoculation 
with S. sclerotiorum in B. napus (Bloem et al., 2012). Glucosi- 
nolates have also been shown to be important for the defence 
response in Brassicaceae (Wittstock and Halkier, 2002). Glu- 
cosinolate was probably associated with S. sclerotiorum resis- 
tance in B. napus (Li et al., 1999). And as discussed above, 
BnaC.IGMT5.a, an indole glucosinolate methyltransferase, is a 
candidate gene for S. sclerotiorum resistance (Wu et al., 2013). 
We therefore conclude that the increased metabolism of 
sulphur-containing compounds, such as GSH and glucosinolate, 
plays a key role in plant defence. 

Specific genes in JA signalling and lignin biosynthesis were 
found only in the R genotype, indicating that these two pathways 
have important functions in the defence response. JA, synthe- 
sized from linoleic acid by a series of enzymes (i.e. LOX, AOS, 
AOC, and OPR), is an endogenous hormone involved in plant 
defence against necrotrophic pathogens (Frendo et al., 2013). 
Using microarray analysis, Yang et al. (2007) and Zhao et al. 
(2009) found that genes associated with the JA pathway were 
induced in response to S. sclerotiorum in B. napus. The sequen- 
tial activation of salicylic acid (SA) signalling and JA signalling has 
also been shown to be important in the defence response against 
S. sclerotiorum in B. napus (Wang etal., 2012b). The 
JA-mediated defence response is thought to be though the 
regulation of protease inhibitor and secondary metabolites, such 
as anthocyanin and glucosinolates (van Dam etal., 2004; 
Mengiste, 2012). 

Lignin in the plant cell wall contributes to disease defence, not 
only by enhancing the mechanical strength of the plant, but also by 
preventing infection and colonization of the pathogen due to its 
hydrophobic properties (Vanholme et al., 2008). Lignification 
protects cotton from the fungal necrotrophic pathogen Verticil- 
lium dahliae (Shi et al., 2012; Xu et al., 2011). In our study, specific 
lignin synthesis-related genes, that is C4H, CCOAOMT, COMT and 
CAD5, were up-regulated in the R genotype, which shows that 
lignin contributes to the response to S. sclerotiorum in B. napus. 
Miedes et al. (2014) summarized and explained the effect of lignin 
on plant pathogen resistance or susceptibility, which was mainly 
through the alternation of cell walls and the production of some 
soluble phenolic compounds. However, a greater understanding 
of the biochemical mechanisms is still necessary to fully understand 
the lignin host-pathogen interactions. 
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The relationship between GSH and other stress response 
genes has been studied in response to abiotic and biotic stress 
(Han et al., 2013; Kumar et al., 2015; Wingate et al., 1988). 
GSH was found to induce biosynthesis of phenylpropanoid, 
lignin, dihydro flavonols pathway, JA signalling, and the 
expression of TFs and various PR proteins (Garg et al., 2013; 
Han etal., 2013; Kumar etal., 2015). These genes and 
pathways were also induced in our study. Therefore, it could 
be speculated that GSH may also play a role in modulating the 
JA and lignin pathways in response to necrotrophic pathogens 
such as S. sclerotiorum in B. napus. GSH may act as a master 
switch in regulating pathogen response in plants. The increased 
levels of GST could enhance GSH activity (Garg et al., 2013), 
and GSTU detected above may play a major role in the response 
to S. sclerotiorum in B. napus. 


Genome-wide prediction 


We performed a genome-wide prediction to evaluate the 
breeding value of a plant based on its genetic loci with rrBLUP. 
rrBLUP performed better in the prediction accuracies for the 
complex traits, such as flowering time, oil content, plant height 
and glucosinolate content in B. napus, than the method of 
BayesB (Wurschum et al., 2014). The predictive ability was related 
to training population, the number of markers and broad 
sense heritability (Owens et al., 2014; Slavov et al., 2014). In 
our study, prediction abilities also increased with the size of the 
training population, but the standard deviation increased also. 
Therefore, it was better to select 60% of the training population 
(145) for genome predictions. When all SNP markers were used, 
the predictive ability was too low (0.27) that the markers failed to 
predict stem resistance. It is well known that GWAS and genomic 
prediction are performed on the same phenotype and same 
genotype data sets, and that GWAS is able to detect markers 
significantly associated with the trait by setting a certain P value, 
such as P < 0.001, P < 0.01 and P < 0.05 (Zhang et al., 2014b). 
When using the significantly associated SNPs with S. sclerotiorum 
resistance produced by GWAS including P < 0.001 (69 SNPs), 
P<0.01 (288) and P<0.05 (1203), the predictive ability 
increased from 0.58 (P < 0.001) to 0.8 (P < 0.05). Therefore, 
these significant markers detected in GWAS could be used to 
predict stem resistance in efforts to breed B. napus lines that are 
resistant to stem rot. Zhang et al. (2014b) also showed that 
GWAS could improve the accuracy of genomic prediction in 
animal breeding, rice breeding and human genetics, when the 
significant SNPs detected by GWAS and located in QTL regions 
are employed in genomic prediction. With the help of genetic 
architecture and population structure revealed by GWAS, the 
efficiency of genomic prediction and further breeding processes is 
increased (Spindel et a/., 2015). 


Experimental procedures 
Plant material and phenotypic evaluation 


A total of 347 B. napus lines were collected from spring, winter 
and semi-winter accessions and cultivated in Southwest University 
of Beibei, Chongqing, China, in 2012 and 2013 (Table S1). Each 
line was planted in two rows of 10 plants per row, with 30 cm 
between rows and a distance of 20 cm between plants within 
each row. 

The S. sclerotiorum isolate was provided by Dr. Jiagin Mei 
from Southwest University and was cultured on potato dextrose 
agar medium (20% potato, 2% dextrose and 1.5% agar) at 
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22 °C in darkness. The detached stem inoculation assay was 
used to evaluate resistance to Sclerotinia at the end of 
flowering time according to Mei etal. (2012), with ‘Zhong- 
shuang 9’ as control, which is known to be resistant to stem 
rot. A 40-cm segment of stem that was 10 cm above the soil 
was cut off the plants. The cut ends were sealed in plastic 
wrap, and the stems were wrapped in moist towels and filter 
papers and placed on a board. Each stem was wounded at two 
points with a 4-mm-diameter punch, and Sclerotinia mycelium 
agar plugs (6 mm in diameter) were placed on each wound. 
The stems were sprayed with water mist, and the board was 
covered with plastic film. The infection temperature was 
maintained at 22 °C. Lesions first appeared 24 h after inocu- 
lation and expanded to the two sides of the stem. After 72 h, 
lesions were watery and encircled the stems. The lesion lengths 
on the stems were recorded 72h after inoculation. The 
resistance assessment (RA) was calculated according to relative 
susceptibility, RA = L/Leontro Where L denotes the lesion length 
of the tested accession and Leontroi represents that of ‘Zhong- 
shuang 9’. The resistance was conducted with two replications 
for each of 347 lines during 2013 and 2014 using the detached 
stem inoculation assay. 


SNP genotyping and in silico mapping of SNPs 


SNP genotyping was performed using the Brassica 60K Illumina 
SNP array (http:/Avww.illumina.com/technology/beadarray-tech- 
nology/infinium-hd-assay.html), and SNP data were analysed 
using Illumina BeadStudio genotyping software. SNP call fre- 
quencies of <0.9 or minor allele frequencies (MAFs) of <0.05 
were excluded in the study. 

To identify the physical position of SNPs in the B. napus 
genome, SNP array sequences were aligned with the genome 
sequences of B. napus (Chalhoub, 2014) using local BLASTn with 
E values of <'1.0E~'°’. SNP positions with a minimum E value and 
a maximum score were selected. SNPs with the same E values and 
scores were removed. 


Relative kinship and haplotype block construction 


The relative kinship matrix of 347 B. napus lines was computed 
using SPAGeDi software (Hardy and Vekemans, 2002). Negative 
kinship values were set to zero. Haplotype blocks were 
constructed using the four gamete rule (Lorenz et al., 2010) 
with Haploview software (Barrett et a/., 2005). The parameters 
were set as follows: P value cut-off of Hardy-Weinberg was 
0.001, MAF was 0.05, and the maximum number of Mendel 
errors was 1. 


Genome-wide association analysis 


The best linear unbiased prediction (BLUP) of a stem resistance 
trait with two replicates over two years was estimated using an R 
script (www.extension.org/pages/61006) based on a linear 
model. An association analysis was performed using the P+K 
model implemented in GenABEL software (Aulchenko et al., 
2007). P values were calculated based on the mm score function 
of GenABEL. The quantile—quantile plot was displayed with -log10 
(P) of each SNP and expected P value. The Manhattan plot was 
displayed using qqman software (Turner, 2014). The threshold of 
GWAS was set to P < 3.23 x 107° (1/total SNPs used). The false 
discovery rate (FDR) was calculated based on the formula: 
[(m x P)/n] x 100%, where m was the number of total SNPs, P 
was the threshold of GWAS, and n was the number of significant 
SNPs. 


Stem resistance prediction 


To evaluate the prediction ability of markers for stem resistance 
among 243 lines (miss rate<10%), prediction accuracy was 
examined using ridge regression best linear unbiased prediction 
rrBLUP (Endelman, 2011) in the R package. The prediction ability 
(r) of the model representing the correlation between the 
predicted and field data was estimated (Slavov et al., 2014). To 
assess the performance of genomic prediction, 20-90% of the 
population was selected at random as the training population and 
markers (P < 0.001, P < 0.01, P < 0.05 and all SNPs) were also 
used. 


Transcriptome sequencing and the identification of 
transcription factor genes 


Five highly resistant (R) and susceptible (S) lines were selected 
from the GWAS population, and two stems were used per line 
and two inoculation sites on each stem. The R lines were WH-57, 
2011-7103, Huayou14, Chuxianbaihua and Youyan2, and the S 
lines were Yangjian8, P685, Guangde138, SWU69 and 07037, 
respectively (F = 48.69, P < 0.0001). Mock-inoculated lines were 
treated with agar plug only. The stem tissues extending 10 mm 
beyond the inoculation point and 1 mm deep were harvested 
48 h postinoculation with a razor blade, and immediately placed 
in liquid nitrogen and stored at —80 °C. Five R and S lines (20 
inoculation points) before and after inoculation were pooled to 
yield four RNA samples, that isRo, Rag, So and Sag. Sequencing 
libraries were generated using the Illumina RNA Library Prep Kit 
following the manufacturer’s recommendations. The library 
preparations were sequenced on an Illumina Hisegq 2000 plat- 
form, and 100-bp paired-end reads were generated. 

The B. napus reference genome was constructed using bowtie 
2 (Broad Institute, Cambridge, MA, USA). Sequencing reads were 
aligned to the B. napus reference genome and then assembled 
using TopHat 2.0.0 and Cufflinks (Trapnell et al., 2012). Gene 
expression levels were estimated using FPKM (fragments per 
kilobase of exon per million mapped fragments). Differentially 
expressed genes (DEGs) between two samples were identified 
with Cuffdiff, based on the criteria FDR<0.01 and |logz (FPKMas/ 
FPKMo)|>2. To identify possible transcription factors (TFs), DEGs in 
the R and S genotypes were aligned to TFs in A. thaliana, which 
were downloaded from the Plant Transcription Factor Database 
(http://planttfdb.cbi.pku.edu.cn/index.php). 


GO and KEGG enrichment analysis of DEGs 


GO enrichment analysis was implemented using BiNGO in 
Cytoscape v3.1.0 (Maere et al., 2005; Smoot et al., 2011), and 
significantly enriched GO terms (FDR < 0.01) were displayed using 
the online tool WEGO (Web Gene Ontology annotations Plot, 
http://wego.genomics.org.cn) (Ye et al., 2006). The enrichment 
of DEGs in the KEGG pathway analysis was conducted in the 
KOBAS2.0 website (http://kobas.cbi.pku.edu.cn/nome.do). To 
analyse the metabolic pathway and functional classification of 
DEGs, expression data were mapped to metabolic pathways using 
MapMan software (Thimm et al., 2004). 


qRT-PCR analysis 


Total cDNA was synthesized from 1 pg RNA according to the 
manufacturer's instructions (Bio-Rad, Foster City, CA). Primers are 
listed in Table S2. PCRs contained 10 uL SSoAdvanced SYBR 
Green Supermix (Bio-Rad), 2.0 uL cDNA, 1 uL primer and distilled 
water to a final volume of 20 uL. The PCR program was as 
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follows: 95 °C for 30 s and 35 cycles of 95 °C for 5 s, followed 
by 56-67 °C (depending on the primers used) for 30 s. Three 
technical replicates were run for test and reference genes. 
Relative expression levels were evaluated using the 2~44¢T 
method. Actin was used as an internal control for normalization, 
and the inoculated group was compared to the mock-inoculated 
control group. 
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